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Abstract 

The (3 — v correlation coefficient, an v , is measured in 21 Na by detecting the time-of- flight of 
the recoil nucleus detected in coincidence with the atomic electrons shaken off in beta decay. 
The sample of 21 Na is confined in a magneto-optic trap. High detection efficiency allows low 
trap density which suppresses the photoassociation of molecular sodium, which can cause a large 
systematic error. Suppressing the fraction of trapped atoms in the excited state using a dark trap 
also reduces the photoassociation process, and data taken with this technique are consistent. The 
main remaining systematic uncertainties come from the measurement of the position and size of 
the atom trap, and the subtraction of background. We find ap v = 0.5502(60), in agreement with 
the Standard Model prediction of ap v = 0.553(2), and disagreeing with a previous measurement 
which was susceptible to an error introduced by the presence of molecular sodium. 

PACS numbers: 23.20.En,23.40.Bw,24.80.+y,13.30.Ce 
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I. INTRODUCTION 



Precise measurements of beta-decay kinematic correlations test the Standard Model of 
the Weak Interaction. The correlation of the momenta of the (3 and the v is sensitive to 
scalar and tensor current interactions which are suggested by some Beyond Standard Model 
theories [l[ 2], [J . Neutral atom optical traps containing short-lived isotopes are an appealing 
source for the study of (3 decay correlations, and several experiments have been performed 
and proposed in such systems 

anaa 

. Recoiling daughter nuclei (energy ~100eV) and 
13 particles emerge from the trap volume with negligible scattering and propagate in ultra- 
high vacuum. Decays occur essentially at rest (temperature <250 /xK) and are localized 
in a volume smaller than 1 mm 3 . The trapped atoms are isotopically pure, and very little 
background activity is present when the detection region is restricted to a small zone around 
the trap. The source location and spatial distribution can be monitored by observing atomic 
fluorescence. The polarization state of the nuclei can be manipulated using optical pumping 
techniques. 

A previous measurement of the beta-neutrino correlation in laser-trapped 21 Na found 
a value 3.6<r smaller than that calculated from the "V minus A" current coupling of the 
Standard Model In that work, the results suggested a dependence on the number 
of atoms: at lower trap populations, the measured correlation coefficient seemed larger, 
and if extrapolated to zero trap population, the result agreed with the calculated value. 
This tentative conclusion is addressed here: repeating the experiment using a more efficient 
detection system, we find convincing evidence that the previous measurement of ap v was 
distorted by events originating from cold, trapped molecules of 21 Na2- A re- measurement of 
ap u has been made with particular attention to minimizing the effect of molecular sodium. 

II. THE BETA-NEUTRINO CORRELATION IN THE ELECTROWEAK STAN- 
DARD MODEL 



In a source with no net nuclear po 
can be inferred from the (3 decay rate 
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by detecting the (3 and recoiling nucleus in coincidence. Here (E e ,p e ) and (E u ,p v ) are the 
(3 and z/ 4-momentum, Eq is the /9 decay endpoint energy, m e is the electron mass, and 
ermi function from [ll|]. The energy distribution of the recoil nuclei is 



F(Z, E e ) is the 

calculable (see [3, 12j) but gives no greater insight into the decay rate dependence, and the 
calculation is most simply performed in terms of integrals over (3 energy. In the allowed 
approximation, fx, the (3 — v correlation coefficient (a^„) and Fierz interference term (&Fierz) 
are independent of E e . Their values are determined by fundamental weak coupling constants 
Ci and G[ {i = scalar (S), vector (V), tensor (T), and axial- vector (A)}, and by the Fermi 
(Gamow- Teller) nuclear matrix elements, Mp (Mqt) |9(. For a mixed Fermi/Gamow- Teller 
beta decay, 

ap v = (|M*| 2 (|<7vr| 2 +|<^| a H<7*| a H^^ (2) 
where 

i = \M F \\\C V \ 2 + \C' V \ 2 + \C S \ 2 + \C' S \ 2 ) + \M GT \ 2 {\C A \ 2 + \C' A \ 2 + \C T \ 2 + \C' T \ 2 ), (3) 
and 



6 Fierz = ±2 y/1 - (ZaY Re [\M F \ 2 {C S C* V + C' S C%) + \M GT \ 2 (C T C* A + C' T C'X)\ C ■ (4) 

In the Standard Model (SM), Cy and Ca are nearly purely real (and we neglect any imag- 
inary component for the rest of our analysis), Cy = C' v = 1, and Ca = C' A ~ —1.27 (from 
experiments), and all other coupling constants are zero. Experimental limits on scalar 
and tensor couplings predicted by some SM extensions are model dependent (depending 
on assumptions about theparity and time- reversal properties of new interactions), and not 
stringent in some cases If present, these couplings would alter ap u and could cause 
non-zero Fierz interference terms. 

Several corrections alter the allowed approximation prediction of a^ v and give E e depen- 
dence to fi, ap v , and ^Fierz at the 1% level. The ground state 21 Na(3/2 + ) decays via two main 
branches to 21 Ne(3/2 + , OkeV) and 21 Ne(5/2 + , 350.7keV). Electron capture contributes a 
0.1% branching ratio. From isospin symmetry of the mirror decay 21 Na (3/2 + )-^ 21 Ne (3/2 + ), 
M F = 1. The Gamow- Teller contribution, A = \CaM gt \, can be determined by the ratio of 
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the ft value to the (0 + — > + ) superallowed transitions. The ft value for the ground state 
transition (3/2 + ) is determined by the half-life ti/2 = 22.49(4) s, and the total decay energy 
3547.07(21) from [l3]. The recoil order corrections were evaluated in Ref. [l^] and include 
order-a radiative, and isospin symmetry breaking corrections. We calculate ft = 4036.7(9.0) 
and _ag u = 0.553(2) for 21 Na decays to the ground-state, which is substantially the same as 
in 14|, but accounts for the most recent decay energy measurement. The decay to the 
(5/2 + , 350.7 keV) state is pure Gamow- Teller, and so ap u (5/2 + ) = —1/3 (from Eqs. [2] and 
[3] when Mp = 0) to a good approximation. The omission of recoil order and order-a ra- 
diative corrections, the excited-state lifetime, and (3 — 7 correlation is inconsequential for 
the decays to the (5/2 + , 350.7 keV) state. The decay branching ratio to this state has been 
re-measured to be 4.74(4)% [ijj], and 4.85(12)% 15], compared to the previously accepted 
value of 5.02(13)%. The new measurements confirm that the branching ratio was not a 
source of error in our previous measurement of anv. 

n n 

The magnetic moments of the 21 Na [16J and 21 Ne 17J] nuclei determine the weak mag- 



netism recoil order term, characterized by &wm = 82.6, assuming the conserved vector cur- 
rent hypothesis. The absence of second-class currents implies that the induced tensor term 
is zero for an isodoublet transition. Order-a radiative effects are included in our calculations 
of the decay phase space according to the prescription in Ref. 18j, which accounts for hard 
bremsstrahlung emission. 



III. (3 - v MEASUREMENT TECHNIQUE 



Much of the experiment is similar to the apparatus and technique described previously in 



20] . Sodium-21 is produced via the 24 Mg(p,a) 21 Na reaction at the 88" Cyclotron at 
LBNL. A beam current of between one and two microamps of 25 MeV protons is incident on 
a natural abundance magnesium oxide target enclosed in a sealed, heated, ceramic oven. The 
21 Na diffuses from thin disks of the target material (ten disks of ~ 40 mg/cm 2 ), and emerges 
from the oven (at about 1000 C) as a neutral atomic beam. The beam is collimated with 
20 mm long by 0.6 mm diameter tubes and with a two-dimensional optical molasses region 
just downstream from the oven. The atoms are slowed using a tapered magnetic field and 
counter-propagating laser beam (Zeeman slower). The slowed fraction of the atomic beam is 
loaded into a magneto-optic trap formed by three perpendicular, retroreflected laser beams 
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of about 10 mW/cm 2 each. The magneto-optic trap (MOT) is located near the center of 
the counting chamber containing two microchannel plates (MCPs) and several electrodes to 
form a focusing electric field of about 1.6 kV/cm in the region of the trap. The microchannel 
plates are each about 87 mm away from the trap along a single axis (x). A sample of 21 Na 
is accumulated in the magneto-optic trap. The f3 decay leaves 21 Ne in a variety of charge 
states through shakeoff and Auger processes 19]. The electric field accelerates the ionized 
recoil nuclei to one MCP biased at about -7 kV (in the +x direction), and the low-energy 
electrons shaken off by the 21 Ne towards the second MCP biased at about +5.6 kV (—x 
direction). A trigger from the electron MCP starts a 3/is coincidence window on a time to 
amplitude converter (TAC), with a stop signal provided by the 21 Ne ions detected by the ion 
MCP. The TAC signal is digitized (along with pulse height and rate information) to count 
the time-of-flight (TOF) spectrum for the recoil ions. The (3 — v correlation can be inferred 
from the TOF spectrum since aligned lepton momenta (caused by the ap v term in Eq. [1]) 
result in larger nuclear recoil energies. An example of the TOF spectrum data is shown in 
Fig. [1], representing about seven hours of data acquisition with traps containing an average 
of 5 x 10 4 atoms. Events from neutral 21 Ne are visible as a peak at about 2000 ns. The 
detection efficiency for neutral recoils is low, since these recoils are unfocussed by the electric 
field, and because of the low intrinsic efficiency of the MCP (which has no cesium iodide 
surface coating as in our previous configuration). Also visible is a peak at 786 ns caused 
by autoionized 21 Na2~ formed in the MOT (see Section IVl below) . During data acquisition, 
a CCD camera acquires images of the MOT approximately every ten minutes. The trap 
population is inferred by measuring the intensity of the fluorescence of the trapped atoms 
using a calibrated photomultiplier tube in current mode. 



IV. GENERATING FIT TEMPLATES 



To interpret the time-of-flight spectra as a measurement of a^, we use Monte-Carlo 
simulations to generate template TOF curves for the two kinematic terms in Eq. [1] (ignoring 
^Fierz) with a>p u = and ap u = 0.553. The recoil ion TOF data are then fit to a linear 
combination of these two template spectra to determine the magnitude of ap v . We construct 
recoil ion momentum spectra by generating random lepton momenta which are consistent 
with the phase-space distribution (Eq. CD), and then obtain the recoil momenta by momentum 
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conservation. The recoil ions are assigned a random initial position consistent with the 
spatial distribution of the trapped atom cloud as determined from camera images of the 
trap. The recoil ion trajectories are calculated using an ion optics program which calculates 
the electric field (E (af, y, z)^j given the applied electrode voltages and a full three-dimensional 
description of their geometry (Simlon version 8.0). A Monte-Carlo event is accepted after 
applying position- and energy- dependent efficiency functions for the microchannel plate 
detector. The Monte-Carlo generated TOF spectra for beta decays to singly charged 21 Ne, 
calculated with ap v = and ap v = 0.553 are shown in Fig. [2j We generate events for 
decays both to the 21 Ne ground state and to the 21 Ne(5/2 + , 350.7 kev) state, weighted by 
the branching ratio. Electron capture events are included with a 0.1% branching ratio. The 
effect of electron capture events in the time-of-flight templates is not resolved in our data 
analysis to measure ap v . 

A. Trap cloud parameters 

Measuring the size and location of the trapped atom cloud with respect to the MCP 
detectors is crucial for determining ap u . For each data run, the position and size of the 
magneto-optic trap were determined from camera images and measured time-of-flight distri- 
butions. The scale of the camera image was determined using calibrating rulers at variable 
distance from the camera lens. The MOT position was measured with respect to locating 
cross hairs (20 micron wires) placed at the entrance and exit optical ports along the viewing 
axis of the camera. The cross-hair camera pixel locations were observed at the beginning 
and end of each data group to correct for any camera movement. This locates the MOT 
in a plane perpendicular to the camera viewing axis (the y' — z plane, where y' is oriented 
at 45 degrees with respect to vertical (+y and +x), with an uncertainty of 0.04 mm. The 
location of the cross-hair based reference point with respect to the microchannel plate was 
determined in a calibration using molecular ions' time-of-flight. The MOT location was var- 
ied along the y and z axes using shimming magnetic fields. Using an electric field a factor of 
two lower than is used to detect (3 decay recoils (to avoid strong focusing of the low-energy 
dimer ions), we measured the TOF of autoionized dimer ions ( 23 Na2~) as a function of trap 
position y and z, relative to the cross hair reference point. For the nearly uniform electric 
field directed along x towards the ion MCP, the TOF has a two-dimensional parabolic de- 
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pendence on the trap's plane position (y, z). The minimum TOF occurs nearly at the center 
of the ion MCP in the y — z plane (verified by ion trajectory computation). The dimer 
ions' time of flight at this location depends only on the x position, i.e. the distance to the 
microchannel plate. The absolute time of flight was determined from the TAC value with 
an off-line calibration of the TAC using a precision digital delay generator. The T = offset 
of the TAC was determined from coincident events in which secondary electrons from the 
surface of the ion MCP, propagated back to the electron MCP, with a delay from the electron 
flight time calculated at AT = —8.9(2) ns. The x position of the MOT was obtained using 
the calculated dimer ions' trajectory. The trap location in the yz plane was measured with 
a position uncertainty of 0.2 mm. The absolute TOF of the ionized dimers at the location 
of the minimum TOF determined the x position to an uncertainty of 0.01 mm, (accounting 
for electric field uncertainty). The distance from the MOT to the ion MCP was verified by 
measuring the TOF of 21 Na^ during a beta-decay data run. In this case, the TOF of the 
radioactive dimer ions was sensitive mainly to the x position, when the y, z position was 
known from the camera image. This measurement gave an x axis position of 86.84(1) mm 
from the ion MCP surface for this run. The x axis distance was also measured using /3-decay 
coincidence events of neutral recoils, 21 Ne°, finding 87.0(3) mm. In this case, the rising edge 
of the TOF spectrum depends on the MOT/MCP distance, independent of the electric field, 
and nearly independent of the value of ap u and the size of the MOT cloud. This technique 
was used only in some of the data runs with resolved 21 Ne° events. 

The spatial distribution of the atoms in the MOT was assumed to be gaussian in the 
Monte-Carlo procedure. The full-width at half-maximum (FWHM) of the MOT spatial 
distribution was determined by a least-squares fit to the camera images. A typical FWHM 
was 0.6 mm, varying by about 30% between different data groups, depending on the MOT 
laser and magnetic field parameters. The gaussian FWHM was up to 30% different along 
different axes of the trap, and a two-dimensional Gaussian fit with an asymmetry parameter 
was typcially performed. During data runs, the MOT position remained constant to within 
0.2 mm, determined by the centroids of the fits to the camera images. To analyze each 
data group, the position (x, y, z) and FWHM of the trap were inserted into the Monte-Carlo 
simulations, and the measured potentials applied to the electrodes and MCPs were used to 
calculate the electric field E (x, y, z) to propagate simulated ion trajectories. 
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V. MOLECULAR SODIUM 



As the MOT operates, molecular sodium (Na^) is generated via photoassociation during 



collisions between cold trapped atoms 



211 ] . Cold, ground state molecule formation, using 



only MOT laser frequencies, has been measured in cesium [22| , while a dedicated experiment 
observed trapping of triplet ground state Cs2 by the quadrupole MOT magnetic field gradient 



acting on the molecular magnetic moment 



231 ] . In sodium, the second molecular excited 



state manifold (asymptotically equivalent at large interatomic separation to the state with 
both partners in the 3 2 P^/2 state) is autoionizing. This creates a low-energy electron and 
a Na^ ion with very low initial momentum. We detect molecular sodium as coincidence 
events in the MCP pair with a strongly peaked TOF corresponding to singly ionized species 
with twice the mass of atomic sodium, originating from the trap location. In our previous 
measurement, Ref. |sj], the result for ap u depended on the number of atoms held in the MOT, 
and this dependence was likely caused by counting beta decay events in which the detected 
recoil nucleus scattered from a molecular partner. 

For beta decay occuring in 21 Na2, the recoil nucleus is created near the molecular partner. 
The scattering potential for 21 Ne- 21 Na is not known, and it would depend on the charge 
state of the 21 Ne recoil. A simple model of the scattering of the 21 Ne recoil nuclei by 21 Na 
suggests a very strong perturbation of the measured value of a^ v . Scattering of the recoil 
nuclear momentum will tend to randomize the momentum direction, and the momentum 
would be shared with the molecular partner, distorting the initial momentum distribution 
of the j3 decay recoil nuclei. The scattering process was modeled using a procedure similar 



to Ref. 



24j . assuming a simple Lennard- Jones interatomic potential for the recoil Ne 



with 21 Na. Monte-Carlo generated recoil momenta were generated at a 23 Na- 23 Na binding 



distance r = 0.28 nm 25j, randomly directed with respect to the scattering potential, 



and then propagated in the scattering potential to determine final momenta. These were 
propagated in the ion optics program to create a TOF spectrum for scattered recoil nuclei. 
The simulated, perturbed recoil TOF spectra were fit with simulated, unscattered template 
TOF spectra to determine the approximate size of the error introduced in a measurement 
of ap u . For both the experimental configuration of Ref. 8J (in which the (3 was detected) 
and the current configuration (triggered by shakeoff electrons) the value of a$ v is reduced 
to 0.6 times the nominal value used in the simulations. The reduction factor depends on 
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the simulated Lennard 
- 23 Na parameters in 



ones parameters, but for a range of realistic values (using the 23 Na 
25]), ap v is reduced by about 0.4 to 0.7. 
The autoionization process in our MOT sample is caused primarily by absorption of 
laser photons to form singly-excited, short-range molecular states, which then absorb a sec- 

lis process is known as 



261 ] to be the dominant 



ond laser photon to excite to short-range, autoionizing states. T 
photoassociative photoionization (PAPI), and was shown in Ref. 
ionization process in experiments using a single laser frequency tuned below the atomic 
resonance. The formation of cold, ground state molecules requires the same intermediate 
step - the population of short-range singly excited molecular states, which then decay via 
spontaneous or stimulated emission to bound ground states. The observation of autoionized 
dimer molecules implies a population of cold, ground state molecular sodium. Cold, ground 
state sodium molecules with a net magnetic moment (e.g. in the 3 £+ state) can be confined 
in the magnetic trap formed by the MOT's magnetic field gradient of 10 G/cm, implying 
a cold, trapped, molecular population. It is difficult to determine the absolute fraction of 
trapped 21 Na which is bound as 21 Na2 by measuring the rate of observed ionized dimers. 
In principle, this would help to determine the size of the systematic error introduced by 
recoil scattering. But this would require knowing both the spontaneous emission and au- 
toionization probabilities from both ground and singly excited molecular states. It is also 
difficult to know the distribution of molecular states involved in the photoassociation and 
autoionization, since the photoassociation occurs with only the MOT lasers, rather than a 
separate, frequency-tunable laser which probes the molecular spectroscopy. We have mea- 
sured the rate of autoionized 23 Na2 dimers as a function of the trapped atomic population: 
this measures the relative cold molecular population of the MOT as a function of atomic 
population, since the formation rate of autoionized dimers must be related to the formation 
rate of cold, ground state, trapped dimers. We find a strong dependence of the dimer ion 
rate per trapped atom (the relevant scaling parameter for perturbing a measurement of ap v ) 
on the population of the MOT, shown in Fig. [31 

VI. DATA ANALYSIS 

The total data set consists of 3.6 x 10 6 coincident electron-ion events taken during eight 
separate runs with a total data acquisition time of about 66 hours. Steady-state trap pop- 



9 



ulations during these runs ranged from 5000 to 83000 atoms, and the statistically weighted 
mean number of atoms in the trap for the entire data set was 5 x 10 4 . Each run was analyzed 
separately, generating new Monte-Carlo templates with different input parameters for trap 
size and location, and different electrode potentials. To analyze the data, we apply the TAC 
time calibration (which corrects for a roughly 2 x 10~ 4 differential non-linearity of the TAC 
and ADC conversion) and add the 4 ns shakeoff electron flight time (with a time spread of 0.2 
ns caused mainly by the size of the cloud) to the raw TAC data to generate a true ion TOF 
spectrum. The amplitude of the background in the TOF spectrum is parametrized using 
measured background data (taken with one MOT beam blocked to prevent the formation of 
the trapped population, but still allowing a slowed radioactive atomic beam into the detec- 
tion volume). Since the TOF is determined from the first ion MCP trigger to arrive after 
the electron MCP trigger, lost true and accidental coincidences must be corrected using the 
average ion and electron MCP trigger rates as described in 27J. This requires a slight cor- 
rection to the TOF spectrum, amounting to roughly 0.3% of the counts in the 21 Ne +1 peak. 
The amplitude of the corrected background spectrum is determined by fitting the trap-on 
beta decay TOF spectrum in regions away from the recoil ion and molecular peaks. The 
amplitude of the background term is then fixed in the fits to the recoil-ion TOF spectrum. 
Another issue in subtracting a background TOF spectrum is the contribution from events 
caused by the coincident detection of a f3 + (or 7) with the electron MCP and a positively 
charged recoil 21 Ne ion. Such events were detected by biasing the electron MCP and its 
collimator electrode to reject shakeoff electrons with less than 500 eV originating from the 
MOT location. We estimate the intrinsic detection efficiency for f3 + to be roughly 3% of 
the efficiency for the shakeoff electrons (accelerated to about 1.3 keV) normally used as the 
trigger. The f3 + (or 7) triggered events select a direction for the (3 + momentum, favoring 
recoil ions with a preferred momentum direction towards the ion MCP, which would give a 
distorted measurement of a^ v in the shakeoff-triggered TOF spectrum. It is difficult to ac- 
curately calculate the contribution to the measured TOF spectrum from these /5 + -triggered 
events, since the detection efficiency of the electron MCP as a function of (3 + energy and 
number of incident shakeoff electrons is not precisely known. The /5 + -triggered events should 
have a slightly enhanced coincidence detection probability since the electron MCP is struck 
by both a f3 + and the shakeoff electrons given off in the decay. We correct the measured 
/3 + -triggered spectrum for acquisition time, the slightly different accelerating potential for 
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the ions, and the estimated trigger efficiency enhancement and subtract it from the recoil 
ion coincidence TOF spectrum before fitting. Completely neglecting the subtraction of the 
/3-triggered component to the TOF spectrum changes the fitted value of ap v by only 0.5%. 

In fitting the data by chi-squared minimization, the following parameters are variable: 
the amplitude of each charge state (+1 to +4), an overall time axis shift (T = offset), and 
the ratio of the two Monte-Carlo template curves which account for the isotropic term and 
beta-neurtino momentum correlation term in the beta-decay phase space. The charge states 
branching ratios are known in principle (measured in [191]). but since the absolute detection 
efficiency is unknown for the electron MCP detector as a function of the number of shakeoff 
electrons impinging on it, the strength of each charge state is effectively arbitrary. Data and 
fit for one ten-hour data run are shown in Fig. HI The residuals have been normalized by 
the statistical uncertainty for each point and binned. 

The x location of the MOT was verified by fitting the data using several template TOF 
spectra, generated for a range of x positions. The minimum \ 2 occurred with a MOT-MCP 
distance D = 86.6(2) mm in Monte-Carlo simulation (compared to D = 87.0(3) (based on 
the 21 Ne° TOF) and D = 86.84(1) (based on the dimer ions and optical reference). 



A. Internal Conversion 



A small correction is necessary to the measured value of ap u because internal conversion 
of the excited-state (5/2 + , 350.7 kev) causes the excited state contribution in each charge 
state of 21 Ne to deviate from the j3 decay branching ratio. Internal conversion results in an 
inner shell vacancy subsequent electron loss, causing higher charge states for excited state 
decays, or effectively larger decay branching ratio to the higher charge states in the data. 
Using the same process described in jg], this gives a correction Aap u = +0.0003(1). 



B. Ionization Dependences 

Since the j3 — v correlation is measured only for daughter 21 Ne that have lost >2 electrons, 
the ionization process could lead to systematic effects. As discussed in Ref. [19(, the ratios 
21 Ne +2 : 21 Ne + and 21 Ne +3 : 21 Ne + show no indication the (5 + or recoil-ion energy influences 
ionization. A rough calculation indicates that nuclear recoil should increase ionization for 
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the fastest recoils by 0.70(35)% [19J and we apply a correction of —0.0033(17) to ap u . 
C. Polarization and Alignment 

The nuclear polarization and alignment of the trapped sodium nuclei have been deter- 
mined to be < 0.2% for very broad ranges of trapping parameters by optical rotation spec- 
troscopy on trapped 23 Na atoms 8j. No correction for contributions from polarized nucleus 
correlations are necessary, although a small uncertainty is applied to the final result for the 
uncertainty in the measurement of the polarization state of the atoms. 



D. Molecular sodium 

To address the problem caused by recoil scattering from molecular partners, data were 
acquired with a range of trap populations. Figure [3] suggests that there should be a depen- 
dence of ap u with trap population, as the molecular population fraction per trapped atom 
changes. With a low trap population, according to Fig. [3j we expect a much lower fraction 
of detected beta-decay events will be from molecular 21 Na2. Data were also acquired using a 
dark MOT technique, in which the usual trap repumping laser is not incident on the trapped 
atoms, except at the edges of the trapping region 28|. This optically pumps the atoms to 
F g — 1 where they no longer interact with the trapping lasers. The average 3 2 P^/2 excited 
state population is greatly reduced, inhibiting the photoassociation process (which requires 
collisions involving excited state atoms). With the dark MOT technique, the detected rate 
of autoionized dimers is reduced by a factor of about 700 relative to the bright MOT. The 
measured value of ap v for dark MOT data (expressed as a fraction of the "Standard Model" 
value asM = 0.553) is a = 1.004(17). This value agrees with a = 0.994(7), the mean of the 
data measured with bright MOTs. The dark MOT data was taken with a trap population 
of 50000 atoms. In Fig. [5], we show the current data set and the averaged data from [8] . The 
perturbation on ap u from molecular recoil scattering can be estimated by fitting these data 
to a curve derived from the measured rate of molecular ions (per atom) as a function of trap 
population, Fig. [3J This fit gives a negligible correction (0.05%) to the current data set for 
an extrapolation to zero trap population; while also supporting the plausibility of such an 
extrapolation for the data measured in 
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E. Results 



Averaging the data from the present data set (weighted by the statistical uncertainty for 
each run) gives a^ v = 0.5498(38). Several corrections are applied to this value to account for 
instrument effects and biases caused by the shakeoff ionization process. These are summa- 
rized in Table [J along with the significant sources of systematic uncertainty. The efficiency 
of the ion MCP varies as a function of ion impact location. This was measured by examining 
the pulse height spectra of the ion MCP for events from 23 Najj~ as a function of ion impact 
location (y, z) on the MCP. The effective detection diameter of the ion MCP is uncertain in 
this procedure, and the manufacturer's specification of the detector's active diameter is used 
to estimate this uncertainty. The effect of a slight decrease in efficiency at the center of the 
plate (0.45%) was modeled using the Monte-Carlo generating procedure to create distorted 
TOF spectra, and then fitting this simulated data with undistorted simulated data. A sim- 
ilar technique was applied to study the effect of uncertainty in the electric field caused by 
uncertainty in measurements of electrode potentials and locations in the ion optics calcula- 
tion. The uncertainties in Tabled] are summed in quadrature, while the corrections are added 
linearly. The result is ap v = 0.5502 (38) (46), where the first uncertainty is statistical and the 
second systematic. Again, this is to be compared with the calculated value ap u = 0.553(2), 
which assumes Standard Model current couplings, and which has an uncertainty limited by 
measurements of the decay parameters (half- life, Q value, and branching ratios). 



F. Conclusions 



In Fig.[6]we plot precise measurements of the f3— v correlation for different nuclear systems 
as a function of the Fermi fraction of the transition, F = (l + \C ^Mar / CyMp | 2 ) . We 
show the residual to the theoretical dependence (the measured ap v minus the "V minus A" 
calculated value ogj,(Std. Model)). A plot of a$ v from Eq. [2] versus F would be a straight 
line with offset —1/3 and slope +4/3. Our result is seen to have precision comparable to 
the measurements in 6 He, 38m K, 32 Ar, and the neutron. In these experiments, the beta 
neutrino correlation is measured essentially from the recoil energy spectrum, meaning that 
the observable is 



s = i ■ , isr , ( 5 ) 



1 + &Fierz- 
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with a and &Fierz as in Eqs. [2] and HI and (Ep) the mean energy of the (3 ± . Combining 
the results of several precise measurements of the beta-neutrino correlation a produces a 
limit on the existence of non-Standard Model scalar and tensor currents, shown in Fig. 
To generate this exclusion plot, we assume "normal helicity" C$ = C' s and Ct = C' T , and 
we assume that Im{Cs) = Im{C>r) = (i.e. time reversal invariance). Each system, by 
virtue of the different contriubtions of Fermi and Gamow- Teller transition strengths, and 
different mean beta energies for the Fierz term contribution, yields a different sensitivity 
to possible scalar and tensor contributions. The one, two, and three standard deviation 
contours are derived from the combined constraints offered by the different measurements 
of a. For comparison, we also include on the graph (as a vertical grey bar very near C$ = 0) 
the allowed region (at one standard deviation) of scalar coupling constants consistent with 



the analysis of superallowed Fermi decays 



291 ] . In these systems, the agreement of the Tt 



values of many beta decays despite different Q values provides a constraint to the size of the 
Fierz term in the decay phase space suggested in Eq. [TJ 

Looking forward, measurements of the beta-neutrino correlation coefficient using laer- 
trapped atoms could likely be improved to near the 0.1% level of precision, which would 
offer even more potent limits on Beyond Standard Model physics. In both our experiment 
and that in Ref. [5|, the uncertainties were limited by detector response calibration, electric 
field calibrations, and measurement of the trapped atom cloud distribution and location (i.e. 
the source function uncertainty). The uncertainty generated by the momentum dependence 
of the recoil ionization probability could be reduced by a more sophisticated treatment of 
the problem than the semi-empirical calculation given in Ref. 19|. The uncertainty in the 
predicted value for ap v (0.4%, dominated by the uncertainties in the half-life and decay 
branching ratio) would then be a limiting factor in interpreting a result as a constraint on 
Beyond Standard Model couplings. 
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TABLE I: Corrections and systematic uncertainties. 



Source 


Correction (%) 


Uncertainty (%) 


Recoil ionization 


-0.4 


0.2 


Internal conversion 


0.54 


0.13 


Polarization and alignment 





0.1 


j3 triggered background 


-0.5 


0.4 


£mcp 


0.25 


0.15 


Electric field and simulation 





0.27 


MCP diameter uncertainty- 




0.16 


Trap FWHM 




0.5 


Trap position (x) 




0.2 


Trap position (y, z) 




+0.3 


Total 


-0.11 


0.85 
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FIG. 1: Recoil 21 Ne TOF data. The time bin width is 0.43 ns. Charge states up to +4 are resolved, 
and neutral 21 Ne is observed as a peak at about 2000 ns. Molecular 21 Na^~ was also detected during 
this run and is the peak near 800 ns. 

Contract No. DE-AC02-05CH11231f. 



[1] N. Severijns, M. Beck, and O. Naviliat-Cuncic, Rev. Mod. Phys. 78, 991 (2006). 

[2] O. Naviliat-Cuncic, T. A. Girard, J. Deutsch, and N. Severijns, J. Phys. G 17, 919 (1991). 

15 




450 500 550 600 650 700 

Ion Time-of-Flight (ns) 

FIG. 2: Monte Carlo simulation of the time-of-flight spectra for 21 Ne + given ap v = 0.553 (calculated 
value) and a$ v = 0, and the difference between these two spectra. There are roughly 4 x 10 6 counts 
in the spectra. 

[3] F. Gliick, Nucl. Phys. A 628, 493 (1998). 

[4] Z.-T. Lu, C. J. Bowers, S. J. Freedman, B. K. Fujikawa, J. L. Mortara, S.-Q. Shang, K. P. 

Coulter, and L. Young, Phys. Rev. Lett. 72, 3791 (1994). 
[5] A. Gorelov et al., Phys. Rev. Lett. 94, 142501 (2005). 

[6] D. Feldbaum, H. Wang, J. Weinstein, D. Vieira, and X. Zhao, Phys. Rev. A 76, 051402(R) 
(2007). 

[7] H. Wilschut et al., Hyperfine Interactions 174, 97 (2007). 

[8] N. D. Scielzo, S. J. Freedman, B. K. Fujikawa, and P. A. Vetter, Phys. Rev. Lett. 93, 102501 
(2004). 

[9] J. D. Jackson, S. B. Treiman, and H. W. Wyld, Phys. Rev. 106, 517 (1957). 
[10] B. R. Holstein, Rev. Mod. Phys. 46, 789 (1974). 

[11] H. Behrens and J. Jiinecke, Numerical Tables for Beta-Decay and Electron Capture (Springer- 

Verlag, New York, 1969). 
[12] O. Nachtmann, Z. Phys. 215, 505 (1968). 
[13] M. Mukherjee et al., Phys. Rev. Lett. 93, 150801 (2004). 
[14] V. E. Iacob et al., Phys. Rev. C 74, 15501 (2006). 



16 






% 



--e 

o 
e 



-$ 



If 



10" 6 ffimJ 

io 4 




Trap Population 



FIG. 3: Rate of detected 23 Na2~ per trapped atom as a function of the trap atomic population. 
The trap population is measured by the MOT's fluorescence. 



FIG. 4: Recoil Ne spectra fit with the MC simulation, with bin width 0.43 ns. Residuals are 
shown with 4.3 ns bins for clarity. 

[15] K. Jungmann, private communication (2006). 

[16] O. Ames, E. A. Phillips, and S. S. Glickstein, Phys. Rev. 137, B1157 (1965). 

[17] P. Raghavan, At. Nucl. Data Tables 42, 189 (1989). 

[18] F. Gliick, Comp. Phys. Comm. 101, 223 (1997). 

[19] N. D. Scielzo, S. J. Freedman, B. K. Fujikawa, and P. A. Vetter, Phys. Rev. A 68, 022716 





300 



400 500 
Ion Time of Flight (ns) 



600 



700 



(2003). 



17 



Scielzo 2004 
Vetter 2007 

Molecular scattering dependence 




~l 1 1 

100 200 300 
Trap population 



400 500x10 



FIG. 5: (Color) a$ v at different trap populations. Data from 8J are also included. The solid line 
is a fit to the expected dependence of perturbation from a population of trapped 21 Na2, causing 
scattering of the recoil momenta. 




T 



T 



I 

0.4 0.6 0.8 
Fermi Fraction 



1.0 



FIG. 6: Deviations of measurements of ag„ from calculated values plotted as a function of the 



^ermi fraction of the ft decay. From left to right, the systems are 23 Ne 
33], 19 Ne Q], 21 Na (this work), 38m K fl], 32 Ar 



32, 



3j, and 18 Ne 



35]. 



30], 6 He 



311 ]. neutron 



[20] N. D. Scielzo, S. J. Freedman, B. K. Fujikawa, I. Kominis, R. Maruyama, P. A. Vetter, and 

J. R. Vieregg, Nucl. Phys. A 746, 677c (2004). 
[21] F. K. Fatemi, K. M. Jones, P. D. Lett, and E. Tiesinga, Phys. Rev. A 66, 053401 (2002). 
[22] A. Fioretti, D. Comparat, A. Crubellier, O. Dulieu, F. Masnou-Seeuws, and P. Pillet, Phys. 

Rev. Lett. 80, 4402 (1998). 
[23] N. Vanhaecke, W. de Souza Melo, B. L. Tolra, D. Comparat, and P. Pillet, Phys. Rev. Lett. 



18 



0.10 



0.05 



0.00 



-0.05 



-0.10 



neutron 2002 


P 

It 


— ! 

i— .... i .... B 





-0.10 



-0.05 



0.00 

Cg/Cy 



0.05 



0.10 



FIG. 7: (Color) Allowed regions for scalar and tensor coupling constants (assuming normal helicity) 
for precise measurements of ap u . The one, two, and three standard deviation contours are shown at 
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